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We investigate the coherent transport properties of a DNA chain on a substrate which is subjected 
to a uniform electric field perpendicular to the surface. On the basis of the effective tight-binding 
model which simulates charge transport through DNA, the transmission coefficient, Lyapunov expo¬ 
nent, and localization length are numerically calculated by using the transfer-matrix method. It is 
found that an isolated extended state may appear at the Fermi level for a certain gate voltage when 
the interaction strength between the chain and the substrate is position dependent but independent 
of the base-pair sequence, leading to the gate voltage induced Metal-insulator transition (MIT). 

Moreover, conductance and current-voltage characteristics are also calculated. The relationship of 
Lyapunov exponent distribution to the current-voltage characteristics is discussed. Two different 
conduction mechanisms are proposed depending on effectively delocalized states and isolated ex¬ 
tended states, respectively. These results may provide perspectives for experimental work aimed at 
controlling charge transport through DNA-based nanodevices. 


1. INTRODUCTION 


Since Eley and Spivey [I| predicted that DNA chain 
could be conductive, various ingenious experiments have 
been designed to measure the conductance of DNA. Some 
of them assert that DNA is an insulator ii , but others 
show the semiconductor or conductor properties [1-Q, 
and even find the evidence of its superconductor behav¬ 
ior 0- Although controversy remains for these conse¬ 
quences, it is believed that DNA molecule has the poten¬ 
tial for conducting, whose conductivity not only depends 
on its nature, e.g. the sequence, the sugar-phosphate 
backbone, and the helix configuration, but also suffers 
crucial impacts of the environment, e.g. temperature, 
humidity, solvent, ions, substrate and external field [^,0. 
In 2001, K. H. Yoo et.al have already found in their exper¬ 
iment that the current though DNA can be altered by the 
gate voltage, that is, DNA molecule is analogous to the 
field effect transistor (FET) M- This phenomenon has 
reappeared several times in the later experiments m- 
M- Theoretically, A.V. Malyshev attributes this phe¬ 
nomenon to the impact of the external field on the helix 
structure [l^. In his model, the DNA chain is suspended 
without any substrate. However, DNA chains are usually 
deposited on the surface of some substrate in most of the 
experiments. In fact, several groups have observed the 
strong interaction between the mica/Si02 surface and 
the DNA chain on it, which does have significant impact 
on the charge transport through DNA pH - l^ . So the 
study of disorder and substrate effect in one-dimensional 
DNA chain is still an engaging topic. In fact, with the 
coexistence of the substrate, gate voltage, and disorder, 
a complicated situation may appear in DNA chain. 

In the present paper, we consider a one- 
dimensional(lD) tight-binding model to investigate 
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the combined effect of disorder, the interaction be¬ 
tween DNA chain and substrate, and the gate voltage. 
The disorder is introduced by random distribution of 
poly(dA)-poly(dT) and poly(dG)-poly(dC) in DNA 
chain. An isolated extended state has been founded for 
a certain gate voltage when the interaction strength 
is dependent on the position of chain but independent 
of the base-pair sequence, which leads to the gate 
voltage induced Metal-insulator transition (MIT). The 
dependance of the transmission spectrum and Lyapunov 
exponent (LE) on the energy, the gate voltage and chain 
length are calculated. The associated behavior is studied 
and the relationship of the results to the distribution of 
Lyapunov exponent is also discussed. 

The paper is organized as follows: In the next section 
we describe the basic formalism in our calculations. In 
section 3 we present the main results and discuss the 
related physical implication. The last section is devoted 
to a brief summary of conclusions. 


2. THE BASIC FORMALISM 


We consider a ID DNA chain with interaction between 
chain and substrate as illustrated in Fig. 1. In the tight- 
binding approximation, the Hamiltonian is described as 


U — ^ tyi (c)^C7i-|_i -j-JL.C.)] 

n 

N 

+ '^[{£bn + Vg)blbn - tbniclbn + H.C.)], (1) 

n—1 

where each lattice site represents a poly(dA)-poly(dT) 
base or poly(dG)-poly(dG) base for n G [1, A*]. The op¬ 
erator cl(cn) creates(annihilates) a charge at the nth site 
of the base pair in the DNA chain (I < n < JV), of the 
left electrode (n < 0), of the right electrode (n > JV + 1). 
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The operator h\{bn) creates(annihilates) a charge at the 
nth site of the substrate (1 < n < N). and re¬ 
spectively, denote the site energy and hopping integral 
in DNA chain, is the site energy in the substrate, tbn 
stands for coupling between chain and substrate, and Vg 
is the gate-voltage. 

On each substrate site there is a highly localized orbit 
formed via the interaction of the base pair, the sugar- 
phosphate structure and the substrate surface. So, it 
is reasonable to consider that different base pairs corre¬ 
sponds to different values of coupling strength t^n and 
side site energy that is, for base pair GC, we have 
ebn = ^bGC, hn = tbGC] for base pair AT, we have 
Cfen = CfcAT, tbn = tbAT ■ Obviously, this is a kind of 
longitudinal short-range correlation analogous to a layer 
structure that may produce a band of extended states 


If the base pair GC (or AT) is distributed randomly 
and independently with a probability p (or 1—p), the dis¬ 
order strength W will depend on p and \eAT — eccl- For 
p = 1 or 0 cases, we have a period chain with two energy 
bands, the highest occupied molecular orbital (HOMO) 
and the lowest unoccupied molecular orbital (LUMO). 
For a fixed \€at — f-Gc\ ^ 0 case, the most disordered 
state is given by p = 0.5 and in this case the trans- 
mittivity of the LUMO is almost zero. So the following 
text only focuses on the HOMO with p = 0.5. 

In the site representation, the Schrodinger equation 
= £'|'0) corresponding to Eq. (1) becomes 


(E — tn — ']_1pD,n — l tn'^ D ,n+l tbn'^S,n: 

(E Cbn ^g)'^S,n — tbn'^D,n- 

Here, il)D,n is the amplitude of the wave function at the 
nth site of the chain(the nth base pair site) and ips,n 
is the amplitude of the wave function at the nth site of 
the substrate. After the decimation of V'S.n [H, , the 

corresponding Schrodinger equation of the Hamiltonian 
(1) is 




Here, e(j is the renormalized on site energy, which is de¬ 
fined as 


— En + 


''bn 


E — {tbn + Vg) 


( 4 ) 


The renormalized site energy T includes the influence of 
substrate on conduction state in DNA chain, and it is a 
function of the gate voltage Vg and the incident electron 
energy E [H- 

The Eq. ([3|) can be rewritten as relations between 
coefficient of adjacent pints in the form of the transfer 
matrix 


''>pD,n+l\ ( '>pD,n \ 

, ^D,n ) ~ 


( 5 ) 


with 


{V-E) t„-i 

tn t„ 

1 0 

For a 1-D chain with length L, the coefficients at one end 
are related to the coefficients at the other end with the 
transfer matrices 



The Lyapunov exponents(LE) are the logarithms of 
eigenvalues of the transfer matrix, and can be calculated 
by using the transfer-matrix method, in which the or¬ 
thonormalization procedure is adopted [ 2 ^. 

In this paper, the transmittivity [ 2 ^ and the localiza¬ 
tion length (LL) [2^,[23 which is the inverse of the LE are 
calculated with the transfer matrix (TM). The current is 
computed with the Landauer-Biittiker formula 




/ = -^y T{E){h-fji)dE, (7) 

here fL/niE) = {1 -b exp[{E ± eVdl‘2. - EF)lkBT]}-^ is 
the Fermi distribution on the left/right lead. Ep and 
Vd are the equilibrium Fermi energy and applied voltage, 
respectively. T{E) is the transmission through N base 
pairs. 


3. NUMERICAL RESULTS 

In the numerical calculations we adopt a DNA chain 
with different length N. The lattice structure is il¬ 
lustrated in Fig. I. The values of the parameters 
are chosen as the following [^: e„ is the average of 
the complementary bases energies, cqg = (ec + ec)/2, 
CAT = {eA + eT)l2^. Here, ec = 7.75eU, ec = 8.87eU, 
eA = 8.24eU, ep = 9.14eU. The site energy on the lead 
Em = 5.36 is the work function of platinum. The site en¬ 
ergy CbGC = lO.leU and tbAT = lleU are slightly higher 
than the base pair energy e„. For simplicity, both the 
hopping integral on the DNA and at the contacts are 
fixed at 0.4eU, to minimize the contact scattering [ 2 ^. 
The hopping integral on the leads is 12eU. The surface 
coupling parameters tbGC = 4.4eU and tbAT = 5.0eU 
which are stronger than the usual coupling between base 
pair and sugar_mhosphate structure [la [l8| . The latter 
is about 1.0eU[3,[^ and is neglected here. Fermi energy 
is set to be 5.36eU through out the paper. 

First, we present the results for a DNA chain with 
different length N. Figure 2 shows the rescaled local¬ 
ization length(RLL) A/TV as a function of chain length 
TV with different gate voltage Vg at Fermi energy. Here, 
A represents the localization length. It can be seen, for 
most values of V/, that the A/TV decreases with increasing 
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FIG. 1: Schematic illustration of extending a DNA chain with substrate to a very long chain. The black circles represent the 
DNA chain, white circles denote substrate, and shaded circles stand for the attached leads. Electronic pathways are shown as 
lines. 


N, while, for some given gate voltages ^ 0.57F and 
Vg 2 ~ 3.8817, A/N is slightly increased with N. The RLL 
A/N decays with power rate as increasing N for localized 
state, while it fluctuates throughout for extended state 
1^ . This difference in scaling behavior indicates that 
there may exist delocalized states for given gate-voltages 
in the thermodynamic limit. 

In Fig. 3(a), we plot the Lyapunov exponent as a func¬ 
tion of the energy E with different gate voltage Vg for 
N = 10®. It can be seen, for all the investigated values of 
Vg, that the smallest Lyapunov exponent always lies in a 
narrow region of energy and appears at different energy 
E. Since the transport properties of the system mainly 
depends on the quasi-particle’s states near the Fermi lev¬ 
els, Fig. 3(b) shows Lyapunov exponent as a function of 
the gate voltage Vg at Fermi energy. It can be seen that 
there exist two valley-like structures ai Vg = Vgi ~ 0.57 
and Vg = Vg 2 ~ 3.88, which are exactly the same gate 
voltages that corresponds to the delocalization behav¬ 
iors in Fig. 2. This valley-like structure indicates that 
there may exist a delocalization-localization translation 


around the valleys. Note that only the points at the bot¬ 
tom of the valleys correspond to the smallest Lyapunov 
exponent. Combined the results of Fig. 3(a) and 3(b), 
we could see, at the bottom of the valleys, i.e. Vgi and 
Vg 2 , that the narrow energy region which corresponds to 
the smallest Lyapunov exponent appears around Fermi 
levels. Considering the scaling behavior of the RLL in 
Fig. 2, we could conclude that these states, which cor¬ 
respond to the smallest Lyapunov exponent at Vg = Vgi 
and Vg = Vg 2 in Fig. 3(b), are isolated extended states 
at Fermi energy. 

In order to understand the origin of these isolated ex¬ 
tended states, Eq. ([3]) should be discussed further. Since 
only two kinds of base pairs G-C and A-T are considered 
in present system, there are two kinds of i.e. and 
s'qc- If parameters satisfy the following conditions, 

^'at = ^GC = ( 8 ) 

E€ [e' -2t,e'+ 2t], 


(9) 
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FIG. 2: (Color online.) The rescaled localization length 
(RLL) as a function of chain length N at E = Ef; red 
line for Vg = O.OOR, green line for Vg = 0.30F, blue line 
for Vg = Vgi « 0.57R; pink line for Vg = 2.0017, light blue 
line for Vg = Vg 2 « 3.8817; yellow line for Vg = 4.0917. 

then eq. © becomes 

(e' - E)tj}D,n = tlpD.n-l + #D,n+l 

which corresponds to an Schrodinger equation for a pe¬ 
riodic chain. In other words, the solution of the above 
conditions corresponds an extended state with energy E. 
So Eq. (IH) and (O are a criteria for finding an isolated 
extended states within our model. The origin of this 
extended state comes from the constructive interference 
between DNA chain sites and substrate sites [s^ . Condi¬ 
tion for Eq. dni) ensures E in the propagation band of the 
renormalized period chain with onsite energy e'. Here, e' 
varies with E according to Eq. ([8]). Both Eq. ([8]) and 
Q are independent of GC concentration p and base pair 
arrangements, which indicates that the existence of the 
extended state will not be impacted by specific base pair 
sequences. This mechanism is similar to the cases in refs. 

MM- 

Substitute the renormalized on site energy e' in eq. 
into the above criteria for isolated extended states, Eq. 
dHl) and (in]) turn into the following forms: 

eVg =-| ±/7, (7 > 0). (10) 

E G[e+- 2t, e'+ -f 2t] U [e'_ - 2t, e'_ + 2t]. (11) 

Here, Ae = eat - tGC, a = Ae“^ - tl^^) - {tbAT + 

^bGc), = efcATCfeGcIAe ^{tbAT^bAT ~ ^bGG^bGc) ~ 



FIG. 3: (Color online.) (a) Lyapunov Exponent(LE) as a 
function of energy E with different Vg', curve a, Vg = 017; 
b, Vg = 0.3017; c, Vg = Vgi « 0.5717; d, Vg = 2.0017; e, 
Vg = I7g2 « 3.8817; fj Vg = 4.0917. The open triangle marks 
Ef = 5.36eI7. (b) LE versus Vg at Ef- The DNA length N 
is 10®. 


7 = (3+ 0^/4:, and = eAT + iLT/(““/ 2 ±\/ 7 -ebAT)- 
Therefore, the relation of the isolated extended state en¬ 
ergy E and the gate voltage Vg is two pieces of paralleled 
line segments with slope e and with the terminal vertexes 
depending on site energies and hopping integrals. When 
E = Ep = 5.36el7, Eq. (fTUl) gives two solutions for Vg'. 
Vg « 0.57(= 17gi) and Vg « 3.88(= Vg 2 ), which means, 
for these two cases, that there is an isolated extended 
state at Fermi level Ep. These results are consistent 
with the previous results in Fig. (2) and Fig. (3). 

Figure 4 shows the transmission coefficients as a func¬ 
tion of energy for different Vg and N. It can be seen that 
there is a transmission band for the chain with differ¬ 
ent gate voltages when the chain length is short. With 
the increase of chain length N, some transmission band 
would disappear, i.e. for Vg = 2 case, while others may 
leave a sharp peak around the value of energy of the iso¬ 
lated extended state, i.e. for Vg = 0, Vgi, and 1^2 cases. 
This implies that there are two kinds of mechanisms for 
transmission. The former depends on the effectively de¬ 
localized states [ 3 ^, which forms in a system shorter than 
its localization length, while the latter depends on an iso¬ 
lated extended state. With the increase of system length, 
the effectively delocalized states disappear, however, the 
isolated extended state is still there. In other words, 
with the increase of N, the transmission band decays but 
leaves a sharp peak around the energy of an isolated ex¬ 
tended state for 17g = 0, Vgi, and 17^2 cases, or vanishes 
for Vg = 2 case. 

In Fig. 5, we plot the conductance of chain with differ- 
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ent gate voltage as a function of chain length N. It can be 
seen that the conductance G decreases rapidly with the 
increase of the length N for cases with Vg ^ IgiorV^ 2 - 
This implies that the number of the effectively delocal¬ 
ized states around Fermi energy for such cases decreases 
rapidly with the increase of length of the system. It can 
also be seen that, for Vg = Vgi and Vg = Vg 2 cases, the 
conductance shows a strong oscillation. The oscillation 
might be caused by the coherent scattering of Bloch’s 
wave at the two contacts of the leads [l^. For Vg = 4.085 
case. Go decays exponentially with some oscillation be¬ 
cause the Fermi level is close to the isolated extended 
state. This shows that conductance mainly depends on 
the chain length N and the energy position of the isolate 
extent states. 



E (eV) 






Wl 

0.57 V 






4 4.5 5 5.5 
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FIG. 4: (Color online.) Transmission coefficient T{E) of a 
chain with different Vg as a function of energy. The curve I, 
II, III, IV, V represent for the cases with N = 1000, N = 
5000, N = 10000, N = 50000, and N = 100000, respectively. 
The probability p = 0.5. The results are averaged over 100 
samples. The Fermi level Ef = 5.36eF is labelled by the 
open triangles. 


In order to consider the combined effect of the gate 
voltage and the chain length, in Fig. 6, we plot I — Vd 
characters of DNA with different length N for different 
values of gate voltage Vg. The current / is averaged with 
100 realizations. It can be seen, with the increase of 
the length N, that, for some kind of Vg^ the current / 
decreases and even becomes zero for some certain values 
of gate voltage, while, for the other kind of Vg around 
Vgi and 1 ^ 2 , the current still exists for the regions of Vg 
around Vgi and Vg 2 when N is up to 10®. This shows 
that transport properties of DNA are suppressed by the 
length. 

The difference between these behaviors of currents lies 
in different conduction mechanisms and can be under¬ 
stood by the concept of the effective delocalized statefs^. 
For the short chain with N < A cases, the quasipar¬ 
ticles are in the effectively delocalized states with high 



FIG. 5: (Golor online.) Gonductance G as a function of DNA 
length N with (a) Vg = Vgi, (b) Vg = V 92 , (c) Vg = OeV, 
O.SOeV, 2.00eV, and 4.09eF. The inset Fig. (d) and (e) show 
the details of the oscillation in figure (a) and (b), respectively. 
The oscillation are very dense in figure (a) and (b). The 
results in figure (c) are averaged over 1000 samples. 


transmittivity spread all over the band and conceal the 
impact of the isolated extended state. This situation cor¬ 
responds to the current peak for the case of Vg = 2 in 
Fig. 6. With further increases of N, the effectively de¬ 
localized states eventually vanish, while the impact of 
the isolated extended state emerges. Therefore, two cur¬ 
rent peak bands form respectively around Vgi and Vg 2 
where the isolated extended state lies around the Fermi 
level. However, the two peaks are not exactly at Vgi and 
Vg 2 , because the Landauer-Biittiker formula is asymptot¬ 
ically equal to / = 2eh~^ lEp-eVd/^ ^ow tem¬ 

perature, and the maximum of current does not require 
the transmission peak exactly at Ep, but the maximum 
area under the curve of T{E) in the conducting inter¬ 
val [Ep — eVd/2,Ep + eVd/2]. However, the larger the 
length N is, the more proximity to Igi and Vg 2 the cur¬ 
rent peaks get. In other words, the effective delocalized 
states will vanish with the increase of chain length, while 
the isolated extended states may still exist. 


4. CONCLUSIONS AND DISCUSSION 

We have investigated the coherent transport properties 
of a DNA chain with the interaction between DNA and 
substrate. Using the tight-binding model and transfer 
matrix technique, it is found that there is an isolated ex¬ 
tended state which is robust against the DNA sequence 
rearrangements. Physically, the origin of such an iso¬ 
lated extended state is the interference between DNA 
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N=100000 


V g (V ) 


mesoscale can be achieved by fabricating single molecular 
FET with long disordered DNA chains (e.g., A — DNA) 
adhering to substrate surface. This mechanism is also 
valid for the general nanowire of binary alloy or con¬ 
ductive polymer interacting with substrate surface. Fur¬ 
thermore, since the principle of complementary nucle- 
obases is a natural transverse correlation of short range, 
this mechanism may be generalized to the ladder model 
of DNA, and may therefore impact the spin selective 
transport through double strand DNA as well [15 , which 
needs further study of both theory and experiment. The 
conductance and current-bias characteristics shows that 
Metal-insulator transition (MIT) could emerge by vary¬ 
ing the gate voltage or the chain length. Two different 
conduction mechanisms attributing to effectively delocal¬ 
ized states and isolated extended states, respectively, are 
proposed. These results may provide perspectives for ex¬ 
perimental work aimed at controlling charge transport 
through DNA-based nanodevices. 


FIG. 6: (Color online.) Total current/ of a chain with differ¬ 
ent length N as a function of the gate voltage Vg and applied 
voltage Vd on the two leads with (a) N = 1000, (b) N = 5000, 
and (c) N = 10000; (d) N = 100000. The temperature is 0.1 
K. The results are averaged over 100 samples. Other param¬ 
eters are the same as those in Fig. 2. 

chain sites and substrate sites. Based on this mecha¬ 
nism, the controllable long distance coherent transport in 
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